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Abstract The emergence and magnetic evolution of solar active regions (ARs) 
of /37(5-type, which are known to be highly flare-productive, were studied with 
the SOHO/MDI data in Cycle 23. We selected 31 ARs that can be observed 
from their birth phase, as unbiased samples for our study. From the analysis of 
the magnetic topology (twist and writhe), we obtained the following results, i) 
Emerging j3"f5 ARs can be classified into three topological types as “quasi-/3”, 
“writhed” and “top-to-top”. ii) Among them, the “writhed” and “top-to-top” 
types tend to show high flare activity, iii) As the signs of twist and writhe agree 
with each other in most cases of the “writhed” type (12 cases out of 13), we 
propose a magnetic model in which the emerging flux regions in a (3^5 AR 
are not separated but united as a single structure below the solar surface, iv) 
Almost all the “writhed”-type ARs have downward knotted structures in the 
mid portion of the magnetic flux tube. This, we believe, is the essential property 
of ARs. v) The flare activity of (3^5 ARs is highly correlated not only with 
the sunspot area but also with the magnetic complexity, vi) We suggest that 
there is a possible scaling-law between the flare index and the maximum umbral 
area. 
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1. Introduction 

A sunspot group with umbrae of opposite polarity within a single common 
penumbra is called the (5-sunspot group (|Kunzel[ll960D . It is well known that ma¬ 
jor flares almost always occur in 5-sunspot groups (|Tanakalll975l[^in and LiggettH 
119871 [Sammis et al.ll2OO01) . Hence, investigation of the formation and evolution 
of the 5-configuration is the key in understanding the major flare activities. 

|Zirin and Liggett] (I1987|l classified the formation of 5-spots in three types 
as follows: i) Emergence of a single complex active region (AR) formed below 
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the surface, i.e., the so-called “island (5-spot”, ii) emergence of large satellite 
spots near a larger and older sunspot, and hi) collision of sunspots of opposite 
polarity from different dipoles. Then Hkhsanov et aZ.I (|2004l) . from their morpho¬ 
logical study of (5-spots, classified the magnetic topology of collisional interaction 
between two major emerging flux tubes into three types, namely “vertical (top- 
by-top) collision”, “lateral (side-by-side) collision”, and “frontal (foot-by-foot) 
collision” (see also llkhsanov and Marushinl I2003|l . Although these researches 
show a variety of ways of (5-spot formation, they did not give us a detailed 
physical view of the formation of (5-spots. This is because their classifications 
were phenomenological. 

Succeeding researches (IKurokawalll987lfTanakairi991IILeka et aZ.[[l996l|van Driel-Gesztelyi et al.\ 
I1997lllshii et a/.iri998IILinton et aLlll999IIFan et al.[ll9991 Lopez Fuentes et 

I2000p have found several physical characteristics of active 5-spots, in that they 
have strong magnetic shear structures in their magnetic neutral line area and 
their opposite polarities rotate around each other. These characteristics were 
interpreted as due to the emergence of a kinked flux rope. Ilshii et oD (|1998ll and 
IKurokawa et al.\ (|2002[l proposed that the 5-spot is produced by the emergence 
of a twisted magnetic flux rope whose strong shear produces high flare activities. 

IPoisson et al.\ (|2013|1 suggested that the 5-spot could be due to the emergence of 
a deformed single tube which has a downward convex structure at its middle part 
with kinking. Thus the twist and writhe of a magnetic rope is now considered 
as a key factor for the development of 5-spots and the associated flare activities 
which sometime show sigmoidal morphology ([Bust and Kum^ I1996|l . 

Although a number of case studies revealed the important characteristics of 
5-spot evolution, we have not yet had a conclusive observational and physical 
view of the 5-spot formation. In this article, we tried to obtain more clues on 
this issue by surveying the 5-spots in solar activity cycle 23 paying attention 
to magnetic helicity (twist and writhe). Among all the 5-spots in Cycle 23, we 
selected those which can be studied from their birth on the visible solar disk 
in order to unambiguously study the magnetic connection between component 
sunspots in a group. Further we limited the samples to the regions, as 
they show the highest flare activity (|Sammis et alA I2()()()p . The magnetograms 
and the continuum images taken with the Michelson Doppler Imager (MDI; 

IScherrer et aZ.l [199511 on board the Solar and Heliospheric Observatory (SOHO; 

[Domingo et aL[ I1995p were used to follow the evolution of the ARs. 

Our study shows that the flare-active 5-spots are mainly formed by the emer¬ 
gence of writhed and twisted magnetic tubes, which appear as a quadrupolar 
magnetic configuration on the photosphere. Even the appearance of a more 
complex 5-configuration can be interpreted as a modification to this basic con¬ 
figuration. On the other hand, the flare activity is low when the magnetic tubes 
in the 5-configuration show insignificant writhe or twist. 

In Section 2, we describe the data used in our study and the analysis method. 

In Section 3, we present our detailed analysis on some representative cases and 
summarize our statistical analysis. Finally, we give discussion on our results and 
our conclusions in Sections 4 and 5, respectively. 
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2. Data and Analysis 

2.1. Region Selection 

To select ARs in Cycle 23, we based our study on the following two data 
catalogues: USAF-MWL (also known as USAF-SOON) 0 and USAF/NOAA 
sunspot dataH. The former data base is compilation of daily solar reports of six 
ground-based observatories. The latter mainly consists of data based on space 
observations after the launch of the SOHO spacecraft. When the type of an AR 
is assigned differently in the two catalogues, we choose the more complex one of 
the region as the representative. In this way we selected 200 fi^5 regions which 
are classified so at least once during their lifetime. Four ARs classified as in 
the catalogues were excluded from our dataset as they showed no sign of i5-type 
in our visual check with SOHO/MDI images. Finally we found and selected 31 
candidates which can be studied from their initial emergence on the visible solar 
disk with reference to the Solar-Monitor website 0. 

2.2. Imaging Data 

To follow the evolutions of ARs, we used the longitudinal magnetograms and 
the continuum images taken with SOHO/MDI with a cadence of 96 min and 
a pixel size of 2.0". We utilized the Geostationary Operational Environmental 
Satellite (GOES) X-ray data to estimate the total flare activity of selected ARs. 
Full-disk Fe xii 195A images of the Extreme ultraviolet Imaging Telescope (EIT; 
IDelaboudiniere et aH[T^55|) on board SOHO taken with a cadence of 12 min and 
a pixel size of 2.6" were used to check coronal structures. The Transition Region 
and Coronal Explorer (TRACE; [Handy et aL[[T999p images in I 95 A, I 6 OOA, and 
5000 A taken with approximately 1 min cadence and 0.5" per pixel were also 
used if available. 

2.3. Image Analysis 

2.3.1. Alignment of AR Images 

A time-series of MDI images of a target AR was aligned by the method proposed 
bv lChael(l2001|) and IChae et aLI(l200ip . We converted original 2 arcsec resolution 
images into 1 arcsec ones, and then applied a non-linear mapping to remove the 
differential rotation effect. With these aligned images, we studied the evolution 
of the target AR, and measured the duration of the (5-state by comparing the 
magnetograms with the continuum images. In this study, we set the brightness 
levels of 0.9 for penumbrae and 0.65 for umbrae relative to the intensity of the 
quiet region (see Figured]). 


iftp://ftp.ngdc.noaa.gov/STP/SOLAR.DATA/SUNSPOTJlEGIONS/USAF.MWL 
^http: / / solarscience.msfc.nasa.gov/greenwch.shtml 
®http://www.solarmonitor.org/ 
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Figure 1. A magnetogram (top, gray scale) and a white light image (bottom) of NOAA 
10314. Green and red contours show the intensity levels of 0.9 and 0.65 of the quiet region 
(QR) intensity and delineate QR-penumbra boundaries and penumbra-umbra boundaries, 
respectively (the middle panel). The AR has ^-configuration in the central part. 

2.3.2. Circulation 

When we closely looked at the evolution of an AR, some sunspots exhibited 
prominent rotation around their centers. We estimated the magnitude of rotation 
by measuring the spatially integrated vorticity around the rotating sunspot in 
the following manner. First, we re-aligned the sunspot location in each sequential 
magnetogram to the center of gravity of the magnetic field strength to remove the 
effect of proper motion of the sunspot. Second, we applied the local correlation 
tracking method(LCT; INovember and Simonl I1988P to these images to derive 
the horizontal velocity field (v) around the sunspot. Then, we calculate the area 
integral of the vorticity, namely circulation (C), by 




( 1 ) 


We trimmed away the edge zone of three-pixel width from the vorticity field 
to remove fluctuating and erroneous values before the final estimation of the 
circulation. The procedure of trimming is done by first making a binary image 
from the original vorticity map, and then applying the erosion function ERODE 
in the IDL (Interactive Data Language) software to the binary image which 
results in the shrinkage of island areas in the binary image by three pixels, 
and hnally masking the original vorticity map with the eroded binary image. 
An example of our method for the vorticity field at 17:35UT on 16 January 
2005 is shown in EigurejH Erom left to right, the velocity field overlaid on the 
magnetogram, the vorticity field, and trimmed-away map of the vorticity field 
are shown. Erroneous vorticity values originate at the border between the real 
velocity field and the quiet region where the velocity of magnetic features is zero. 
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Figure 2. Erosion operation applied to the vorticity field at 17:35UT on 16 January 2005. 
Prom left to right are shown a magnetogram, a vorticity field, and the edge-trimmed map 
of the vorticity field, respectively. The velocity field (yellow arrows) and neutral lines (green 
contours) are overlaid on the magnetogram (gray scale). 


The erosion operation can eliminate the contaminations from the real rotational 
elements around the target sunspot effectively as shown in Figure [21 

According to the mathematical definition of circulation, positive (negative) 
circulation corresponds to counter-clockwise (clockwise) rotation. Note that the 
signs are opposite to those of the magnetic twist described in Section [^.3.41 

2.3.3. Tilt Angle of 6-Spots 

The temporal variation of the tilt angle {i.e. the inclination angle of the line con¬ 
necting the bipoles) tells us the orientation of the kink in an emerging magnetic 
tube. Figure [3] schematically shows the case of a tube with left-handed writhe; 
upward-kinked (left panel) and downward-kinked (right panel). 

When they emerge through the photosphere, the tilt angle rotates in the 
clockwise (CW) direction in both cases. The distance between the footpoints 
increases for the upward-kinked case, while it decreases for the downward-kinked 
case. When the tubes with right-handed writhe emerge, the rotation is in the 
counter-clockwise (CCW) direction and the distance between the footpoints 
changes according to the kink direction as in the left-handed case. Thus we 
can infer the orientation of the kink in the emerging magnetic tubes from the 
temporal variations of the inclination angle and the distance between the two 
spots. 


2 . 3 . 4 . Sign of Magnetic Helicity 


A “magnetic tongue” pattern is proposed as a useful parameter to give the sign of 


magnetic helicity of the emerging field ( [Lopez Fuentes et aL[l2()()()l [Lopez Fuentes et a/.[ | 
120081 ILuoni et al.[ 120111) . When a twisted tube emerges above the photosphere, 
the spatial distribution of photospheric magnetic field has an elongated shape 
like a tongue. The tongues of different polarities show anti-parallel orientations 
along the neutral line. Since the orientation of a tongue pattern depends on the 
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Figure 3. Illustrations of the emergence of tubes with left-handed writhe in the cases of 
upward-kinked and downward-kinked configurations. In both cases, the inclination angle of 
the line connecting the two footpoints rotates clockwise and their separation increases or 
decreases according to the orientation of the kink. 

sense of twist of the emerged tube, we can use this parameter as a proxy of 
the sign of helicity (twist) of the emerging magnetic tube as shown in Figure [4] 
(see also Figure 1 of ILuoni et al. 1 I^HTTIl . The circulation of a sunspot is another 
measure of the sign of helicity of the magnetic flux tube, and consistency can be 
checked by comparing of these two signs. 

To avoid confusion, we here reiterate the definitions for the handedness; the 
sign of the left (right)-handed twist or writhe is negative (positive). 


2.3.5. Flare Index 

Several authors have used the soft X-ray flare index (FI) as a proxy of the 
activity level of ARs, as given by {e.g. lAntaloviil 1199611.Toshi and .Toshil 120041 
IAbrainenkoll2005l|Jing et aL[ 1200611 



k 


3 


where me, tum, and mx are the GOES soft X-ray peak intensity magnitudes 
(from 1.0 to 9.9) of C, M and X-class flares. The indices i, j, and k designate 
flares in each class per unit time period. For an AR, two kinds of periods were 
considered in this article; those in which the AR was visible on the disk and 
the AR was in the (5-state. We do not consider small flares less than the C-class 
in this study because they tend to be frequently hidden by major flares or are 
indistinguishable from the background level. 

3. Results 

The (5-spots studied in this article can be classified, from the point of view of the 
emergence morphology, the magnetic connection, the helicity of sunspots, and 
the flare activity, into the following groups: i) top-to-top collision, ii) writhed 
loop, iii) quasi-/?, and iv) the others. The majority of (5-spots in our sample 
are classified into the first three groups. In the next subsection, we select and 
describe the features of a representative example of each group to show our 
classification scheme. 
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Figure 4. Sample images of magnetic tongues. When the head parts of the “Yin-Yang” 
structure forming the magnetic tongues show clockwise or counter-clockwise orientation, they 
indicate negative or positive magnetic helicity, respectively. 


3.1. Case Studies 

3.1.1. Evolution of NOA A10720 

NOAA10720 was the most flare-productive region in our samples. It had the 
longest (5-duration (longer than 213 h), and showed rapid expansion since its 
emergence. 

The magnetic and morphological evolution of the region was as follows. First, 
an emerging flux region (EFRl) emerged as a simple /3 magnetic configuration 
with almost east-west orientation obeying Hale’s law. Then it expanded rapidly 
in size and evolved to the (5-configuration on 12 January 2005 as shown in 
Figure [5] Another new highly sheared EFR (EFR2), which disobeyed Hale’s 
law, appeared on the neutral line as if it tried to penetrate into the opposite 
polarity portion of EFRl at around 16:00UT on 13 January 2005. 

In Figure [SJ we cannot see any signature of writhe in EFRl, as its tilt showed 
no systematic rotational movements. The magnetic tongues of EFRl in the 
initial emerging phase indicated negative magnetic helicity as shown in the first 
row of Figure [S] This means that the tube emerged with a left-handed twist. 
Unfortunately, the tongue pattern could not clearly be seen after the vigorous 
and complex development of EFRl. 

On the other hand, the tilt angle of EFR2 showed a slow CW rotation, 
suggesting the emergence of a loop with left-handed writhe. The twist of the 
magnetic tubes can be estimated by the circulation around the footpoints of the 
magnetic tubes, i.e. sunspots. 

Now let us see the distribution of vorticity around sunspots of p2 and n2 shown 
in Figure ini By specifying the integration area with a circle or an ellipse, which 
covers the target sunspot as compactly as possible, we estimated the circulation 
over the area covering each sunspot. The reliability of the circulation value was 
estimated by changing the size of the calculation area. The area was enlarged 
or reduced by 2 arcsec for this purpose. Besides the procedures described in 
Section 12.3.21 we masked the opposite polarity areas before the integration to 
avoid contamination. The temporal evolution of the circulation of each sunspot 
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is displayed in Figure [T] Both sunspots showed generally positive circulation 
although they showed negative values at times. When we closely looked at the 
vorticity map and the motions of small sunspots in the neighborhood of the 
n2 spot for a day on 16 January, we found that the rapid sliding motions of 
small sunspots relative to the static environment were recorded as false negative 
circulation. If we take this contamination into account, we can assume that 
the circulations around the p2 and n2 spots were mostly positive during the 
evolution, which means that the footpoints of the magnetic tube of EFR2 showed 
a CCW rotation during their evolution. Therefore, EFR2 can be considered as 
an emerging loop of a left-handed twist. 

The timing of the flaring in the region is shown in Figure [5] After the 
emergence of EFR2 on 13 January, many flares have occurred, including three 
X-class flares. High activity triggered by the emergence of EFR2 reminds us 
of the top-to-top collision model proposed by llkhsanov et al.\ (I2()()4|l . The flare 
morphology gave us additional information on the magnetic helicity of EFR2. 
The SOHO/EIT 195 A image in the left panel of Figure IHl shows an inverse-S 
shaped structure over EFR2 after the XI.2 flare that peaked at 0:43UT. The 
inverse-S sigmoid morphology indicates the left-handed twist of the magnetic 
tube. The TRACE 1600 A image in the right panel of Figure IHl shows bright loops 
connecting the two polarities of EFR2. IGno et oI\ (|2013ll argued that negative 
helicity was accumulated in EFR2 by a model of non-linear force-free fields, 
suggesting the emergence of a tube with a left-handed twist. 

Thus, we conclude that the significant flare activity of NOAA 10720 is due 
to the emergence and the intrusion of the twisted and writhed EFR2 into the 
magnetic neutral line of the pre-emerged EFRl. This region is a typical example 
of the top-to-top collision formation of an active <5-type spot. 
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Figure 5. SOHO/MDI magnetograms (left) and white-light images (right) of NOAA 10720 
corresponding to subsequent days of observations. Identified two EFRs are shown in the fourth 
row. The field of view (FOV) of each panel is 300 X 200 arcsec. 
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Figure 6. White light images (left), magnetograms (middle), and vorticity fields (right) of 
the positive (top row) and negative (bottom row) polarity areas of EFR2, respectively. The 
size of each FOV is 100 x 100 arcsec, and its center is at the center of gravity of magnetic 
flux distribution. In the middle panels are overlaid the velocity fields (yellow arrows) and the 
neutral lines (green contours). The areas used to compute the circulation are shown by a circle 
of 28 arcsec radius in the top row and by an ellipse in the bottom row whose major and minor 
axes are 24 and 18 arcsec, respectively. 


SOLA: DELTA_in_Initial_phase_3rd_revised.tex; 24 July 2015; 1:05; p. 10 



























Evolution and Flare Activity of <5-Spots in Cycle 23 


XI.2 X2.6 X3.8 




Figure 7. Time variations of circulation (thin solid lines) and magnetic flux (dashed lines) 
in sunspots p2 (top) and n2 (bottom) in NOAA 10720. The thick solid line shows the running 
average of circulation over five consecutive data points in each panel. The error bars indicate a 
possible range of the circulation values when the size of the analysis area is changes (see text). 
Vertical lines represent the times of the X-class flares. The magnetic flux values are evaluated 
for the full FOV of Figure [6| 
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Figure 8. Time evolution of umbra (dashed), penumbra (dash-dotted), and total sunspot 
(solid) areas in NOAA 10720. Vertical lines indicate X-class (brown solid), M-class (orange 
dashed), and C-class (blue dotted) flares, respectively. The horizontal thick green line at the 
top represents the duration of the (5-state. The triangle at the bottom indicates the central 
meridian passage of the region. 
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Figure 9. Left: A SOHO/EIT 195A image (in negative print) of NOAA 10720 taken at 
01:47 UT, 15 January 2005, after the XI.2 flare peaked at 01:02 UT. White arrows indicate 
two polarities of EFR2. An inverse-S shaped structure is clearly seen. Right: A TRACE lOOOA 
image taken at 03:51 UT on 15 January 2005 between the two X-class flares. A bright structure 
existed between the two polarities of EFR2 even in a relatively quiet period of flare activity. 
Overlaid contours (red for positive and blue for negative polarities) show the line-of-sight 
magnetic fields of 1500 gauss (G) and 100 G (only in the left panel). The FOV of both images 
is 300 X 300 arcsec. 
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3.1.2. Evolution of NOAA10826 

A day after the emergence of 29 November 2005, NOAA10826 formed a 6 - 
configuration and showed high flare activity in the i5-duration of 56 h. 

Its magnetic and morphological evolution is shown in Figure 1101 Two EFRs 
emerged in succession along the east-west line. The ^-configuration was formed 
between the following sunspot of EFRl and the preceding sunspot of EFR2. 

The magnetic tongue patterns of EFRl and EFR2 showed the sign of negative 
helicity (the top row of Figure [TOl), suggesting the emergence of magnetic tubes 
with a left-handed twist. The tilt of the ^-configuration rotated in the CW 
direction, with an average rotational speed of 35° day“^, which is similar to the 
speed of 30° day“^ measured in NOAA10314 bv IPoisson et al.\ (I2013D . As the 
separation between the umbrae of the i5-spot decreased in the first stage, this 
(5-configuration behaved as the emergence of a tube with left-handed writhe and 
downward kink. As the twist had the same sign as the writhe, we suggest that 
the tube was kinked downward due to transformation into the writhe from the 
twist in the layers below the photosphere. 

The magnetic twist of EFRl and EFR2 was studied by the circulation at the 
footpoints of EFRs. We have selected sunspots nl of EFRl and p2 of EFR2 for 
the evaluation of circulation, as they were isolated and free from possible con¬ 
tamination. From the temporal variations of the circulation shown in Figure [TTl 
we conclude that the circulations in both sunspots were predominantly positive, 
which means that the sunspots in the two EFRs rotated in the CCW direction, 
i.e. the tubes emerged with a left-handed twist. 

High flare activity in this region started just after the formation of the 6 - 
conflguration, producing 17 flares in total (4 M-class and 13 C-class) during its 
disk passage and 11 flares (3 M-class and 8 C-class) during the (5-configuration 
(Figure [T2|). The flare morphology gave some insight on the topology near the 
neutral line of the (5-configuration. Figure [T3] shows the case of an Ml.O flare. We 
can see that the bright flare strand had a highly sheared structure with a left- 
handed twist. Since the bright strand was probably produced by the magnetic 
reconnection in the corona, we can infer that the magnetic field near the neutral 
line had left-handed helicity before the flare explosion. 

All the direct or indirect evidence indicates that the EFRs in this region had 
the same left-handed helicity. As the magnetic fluxes of sunspots nl and p2 were 
of the same order of magnitude (Figure [TT]), it is natural to assume that the two 
EFRs were not physically separated but connected like a single writhed structure 
as in the case of NOAA 10314 reported bv IPoisson et al.\ (|2013p . Therefore we 
classify this region as the “downward knotted” type. 
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Figure 10. SOHO/MDI magnetograms (left) and white-light images (right) of NOAA 10826. 
Identified two EFRs are shown in the second row. The FOV of each panel is 300 X 200 arcsec. 
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Figure 11. Time variations of circulation (thin solid lines) and magnetic flux (dashed lines) 
in sunspots nl (top) and p2 (bottom) in NOAA 10826. The format is the same as in Figure [3 
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Figure 12. Time evolution of umbra (dashed), penumbra (dash-dotted), and total sunspot 
(solid) areas in NOAA 10826. The format is the same as in Figure [8] 
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Figure 13. A high-pass filtered TRACE 195 A image (in negative print) of NOAA 10826. 
Contours indicate line-of-sight magnetic fields of 500 G (thick) and 100 G (thin) in positive 
(red) and negative (blue) polarity areas. After the Ml.O flare, highly-sheared bright structures 
appeared on the 5-part of the region. 

3.1.3. Evolution of AR NOAA10050 

NOAA10050 emerged on 26 July 2002. Although the region developed due to 
successive emergence of several EFRs, the flare activity was not high. A few 
(5-configurations formed and disappeared within their lifetime of less than a 
day at different locations in the region. See Figure [T3] for the magnetic and 
morphological evolution of the region. 

The FFRs in the region emerged with their axes in the east-west direction. 
As they evolved, the two opposite polarities of each FFR were separated with 
each other and sunspots of each polarity merged into a dominant center. The 
evolution was very similar to the bipolar emergence model of IZwaanl (|1978|) . 
The 5-configuration formed transiently when a pair of small opposite polarity 
elements incidentally came together. We found no signature of writhe pattern in 
their motion. 

The temporal evolutions of circulation around the dominant polarities are 
shown in Figure [Tsl The preceding negative polarity changed its circulation 
from negative to positive significantly in the early phase of the evolution. After 
several changes in its sign, the positive circulation became dominant in the 
later phase. On the other hand, the following positive polarity always indicated 
negative circulation during the period. The circulations of the two polarities had 
a tendency to show opposite directions. Furthermore, the evolution of circulation 
was influenced strongly by the successive merging of small sunspots. Therefore, 
we may conclude that the actual twist of the dominant sunspots was not high. 
The flare activity in the 5-state of this region was very modest as shown in 
Figure [m 

The evolution of the region is very similar to the /3-type and we classify the 
region as the “quasi-/?” type. 
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Figure 14. SOHO/MDI magnetograms (left) and white-light images (right) of NOAA 10050. 
The FOV of each panel is 350 x 200 arcsec. 
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Figure 15. Time variations of circulation (thin solid lines) and magnetic flux (dashed lines) 
in the preceding negative polality sunspot (top) and the following positive polarity sunspot 
(bottom) in NOAA 10050. The format is the same as in Figure 1111 Dotted vertical lines 
represent C-class flares. 
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Figure 16. Time evolution of umbra (dashed), penumbra (dash-dotted), and total sunspot 
(solid) areas in NOAA 10050. The format is the same as in Figure [8] 

3.2. Statistical Study 

In our case studies of (5-spots, we notice high flare activity in magnetically 
complex (5-spots. In this subsection, we will statistically analyze the dependence 
of the flare activity on the characteristic parameters of (5-spot regions. 

3.2.1. Parameters 

In Table 1, we summarize several flare indices and the parameters which char¬ 
acterize the 31 ARs we studied. The peak soft X-ray flux (XR^ax) and FI are 
the flare indices, while the maximum sunspot area (5'niax), the maximum umbral 
area (5'umax), (5-term (the duration of the (5-state), the handedness of twist and 
writhe, and the type of (5-spot formation are the characteristic parameters of 
ARs. We will comment on the values of the listed parameters in the following 
paragraphs. 

An AR develops and changes its magnetic type and sometimes passes through 
the west limb keeping its high activity. To eliminate possible ambiguity and 
uncertainty in the estimated parameters, we only considered the events that 
occurred in longitudes less than W70°. The values of FI were thus counted 
and listed in Table 1. The values in parentheses are for the flares that occurred 
during the (5-state of the region. The identification of the 5-state was also limited 
to longitudes less than W70°, so that it gives a lower limit in some cases. 

The 5-state in an AR was not always continuous in time. If there appeared 
multiple 5-states intermittently in an AR, we summed all the time spans of 5- 
periods to obtain the “5-duration” (or “5-term”) in this study. When there was 
a gap longer than several hours in the available MDI continuum images, we 
also used TRACE white light images to estimate the 5-duration. However, there 
remained the uncertainties ranging from seven to 31 hours, which is regarded as 
estimation errors for the 5-duration. 

As seen in our case studies, the magnetic field distributions in fi'yS ARs are 
complex and the estimation of the magnetic twist from the circulation is some¬ 
times ambiguous. However, it is well known that the preceding sunspots have 


SOLA: DELTA_in_Initial_phase_3rd_revised.tex; 24 July 2015; 1:05; p. 19 













Takizawa and Kitai 


more rigid structure and longer lifetime than the following sunspots ( |van Driel-Gesztelyi and Petrovay] [ 
119901) . while the following sunspots tend to be affected by turbulent motions in 
the surroundings which cause the fragmentation of the flux tube (|Fan et a,l.[ 

I1993|l . Therefore we mainly measured the circulation of the preceding sunspots 
and adopted their handedness of twist as more reliable representatives of the 
ARs. To decide the sign of writhe, we considered the evolution of the 5-part as 
described in Section B.3.31 

The flare that corresponded to the peak X-ray flux represents only a single 
event of an AR through the observational period and the event did not always 
occur in the 5-state. On the other hand, FI can be measured over the definite 
period of the 5-duration. Therefore, we regarded FI as a better parameter of 
flare activity in the 5-state than XR^ax- 
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Table 1. Summary of parameters for 3lARs studied. The upper part of the table is 
for ARs in the northern hemisphere, while the lower part is for those in the southern 
hemisphere. Here, XRmax means the peak X-ray flux in the observed period of each 
AR. Flare Index indicates the value integrated over the observed period (the value in¬ 
tegrated over the (5-state is shown in parentheses). 5max means the maximum sunspot 
area in the observed period. S'umax means the maximum umbral area in the observed 
period (the value in the 5-state is shown in parentheses). 5-term shows the integrated 
duration of the 5-state. The twist and writhe are indicated by R (right-handed) or L 
(left-handed) . The “Type” column show the emergence type of ARs; TT (top-to-top), 
DK (downward knotted), UK (upward knotted), and QB (quasi-/3), respectively. 


AR 

XRmax 

(Wm-2) 

Flare Index 

‘S'max 

(Mm^) 

‘5'umax 

(Mm^) 

(5- 

term 

(h) 

Twist 

Writhe 

Type 

9165 

59.0 

257.0 

(42.2) 

1851 

291 (280) 

24 

L 

L 

UK 

9511 

120.0 

275.9 

(247.5) 

714 

114 (114) 

20 

R 

R 

DK 

9678 

20.0 

103.1 

(40.3) 

3046 

623 (584) 

26 

L to R 

L 

DK 

9901 

40.0 

53.1 

(53.1) 

2003 

283 (283) 

188 

— 

R 

— 

10412 

9.8 

22.8 

(6.0) 

944 

170 (160) 

27 

R 

R 

DK 

10488 

19.0 

191.8 

(191.8) 

6786 

1789 (1789) 

158 

L 

L 

DK 

10564 

110.0 

219.0 

(209.0) 

3716 

659 (659) 

no 

— 

— 

TT 

10696 

250.0 

1114.9(1062.1) 

3841 

829 (829) 

192 

L 

L 

TT 

10720 

710.0 

2299.3 (2299.3) 

6906 

2147(2147) 

213 

L 

L 

TT 

10956 

2.9 

3.9 

(2.9) 

1054 

151 (151) 

112 

— 

— 

— 

8506 

3.3 

9.1 

(0.0) 

1285 

245 (211) 

8 

R 

— 

QB 

8926 

2.3 

96.6 

(41.6) 

805 

152 (152) 

52 

L 

L 

DK 

9494 

10.0 

59.9 

(26.7) 

1105 

251 (159) 

32 

R 

R 

DK 

9775 

22.0 

115.3 

(11.0) 

1415 

144 (144) 

32 

L 

L 

DK 

9900 

4.4 

11.2 

(0.0) 

633 

75 (67) 

10 

— 

— 

— 

9904 

9.2 

8.2 

(0.0) 

247 

52 (52) 

30 

— 

L 

— 

9906 

26.0 

58.1 

(48.5) 

2907 

540 (540) 

114 

R 

R 

DK 

10017 

150.0 

282.5 

(174.4) 

2298 

341 (341) 

34 

R 

R 

TT 

10050 

8.6 

57.5 

(6.5) 

3408 

627 (600) 

38 

— 

— 

QB 

10119 

8.8 

54.9 

(14.0) 

3192 

541 (541) 

92 

L 

— 

QB 

10137 

40.0 

174.7 

(170.0) 

1330 

252 (252) 

88 

L 

L 

DK 

10226 

68.0 

231.6 

(124.0) 

2812 

578 (578) 

no 

L 

R 

DK 

10314 

150.0 

511.1 

(443.9) 

2367 

635 (635) 

94 

RtoL 

R 

DK 

10417 

5.2 

39.7 

(6.9) 

1405 

199 (158) 

15 

— 

— 

QB 

10456 

4.4 

6.6 

(2.2) 

772 

101 (60) 

14 

— 

— 

QB 

10551 

3.2 

6.9 

(0.0) 

1531 

221 (221) 

38 

— 

— 

QB 

10591 

12.0 

12.0 

(12.0) 

412 

82 (82) 

32 

R 

— 

— 

10798 

56.0 

84.8 

(84.3) 

2931 

461 (461) 

65 

R 

— 

TT 

10826 

78.0 

204.7 

(200.9) 

2152 

390 (390) 

56 

L 

L 

DK 

10848 

41.0 

14.9 

(0.0) 

1466 

181 (181) 

21 

— 

— 

— 

10865 

14.0 

30.0 

(0.0) 

2532 

568 (285) 

38 

— 

— 

QB 


SOLA: DELTA_in_Initial_phase_3rd_revised.tex; 24 July 2015; 1:05; p. 21 







Takizawa and Kitai 


S 8000 

cC 

m 7000 
cb 

pH 6000 

g 5000 

I 4000 

>3 3000 
cS 

g 2000 
g 1000 


1.0 2 2.5 3 3.5 4 




1.6 2 2.5 3 3.6 4 


Figure 17. 5max (left) and FI (right) versus the penumbral-umbral radius ratio in the 5-state. 
3.2.2. Flare Activity and Sunspot Area 

ISammis et al.\ (|2000|1 statistically studied the relation between the flare activity 
and the maximum sunspot area for the ARs in Cycle 22, and found that ARs 
with larger maximum sunspot area produced higher peak X-ray flux. Especially 
they showed that /37i5-ARs had the largest maximum sunspot area and hence 
produced the strongest flares. Now we will check the relation between the FI 
(the proxy of flare activity level) and the maximum umbral area for our sampled 
/?7(5 -ARs in Cycle 23. 

First, we confirm that the relation found bv ISammis et al] (120001) also holds 
in our data. The physical meaning of the relation is that the sunspot area is 
a measure of the total magnetic flux of the region and hence represents the 
magnetic energy content of the region. The ARs with large content of magnetic 
energy can produce large release of energy as flaring. As we do not have the 
magnetic flux values for all the ARs in Cycle 23, we have to use the measured 
area as a proxy parameter for the magnetic flux. The line-of-sight magnetic data 
by SOHO, for example, cannot provide a reliable estimate of the magnetic flux, 
without the correction for the projection effect and the field inclination to the 
vertical. We took the umbral area as a better proxy parameter of the magnetic 
flux than the total sunspot area including the penumbral area, because of two 
reasons. The first is that larger sunspots have smaller contribution of penumbra 
to the total sunspot area (left panel of Figure flTll and contain more magnetic 
flux in the umbral area ((Antaloval 119911 l.lin et aO I2()()fi|l . The other empirical 
reason is that FI does not depend so strongly on the penumbra-umbra radius 
ratio when this ratio is high (right panel of Figure flTll . 

Next, we notice large dispersion in FI in Figure 1181 For a given maximum 
umbral area, there appears a variety of ARs from highly flare productive ones 
to less productive ones. In Figure [THl ARs are roughly split into different 
emergence types, suggesting that the emergence mode is another key parameter 
to characterize the flare productivity. By comparing the two panels in the figure, 
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Suniai in Observation (Mm^) in 6-state (Mm^) 


Figure 18. Flare index versus maximum umbral area in a logarithmic scale. The left panel 
is for the data in the observed period and the right panel is for those in the 5-state. Six 
ARs which showed no flare activity in the 5-state were excluded. Different symbols indicate 
“top-to-top” (red circles), “writhed” (orange squares, both DK and UK types in Table 1), 
“quasi-/?” (blue triangles), and “others” (gray diamonds), respectively. The error-bars were 
assigned by considering the effect of data gaps (see text). The results of non-linear least-squares 
fitting (power-law index /cnl) and classical least-squares fitting (power-law index kc) are shown 
by the solid and the dashed lines, respectively. 


we can see that the separation of the three groups in the right panel is clearer 
than in the left panel. Since the clear separation into groups may give hints of the 
i5-state, we will mainly concentrate to the flare activity and the area variation 
in the (5-state hereafter. 

3.2.3. Flare Activity and Emergence Type 

Figure [m shows the histogram of (Sumax and FI distributions grouped according 
to the emergence type in the (5-state. We combined the groups DK (downward- 
kinked) and UK (upward-kinked) in Table 1 as the “writhed” type. The bin 
sizes of S'umax and FI are I0°-^(~ 2.51) and 10°-^(~ 5.01), respectively. The 
grouping follows the decreasing order of magnetic complexity {i.e. “top-to-top”, 
“writhed”, “quasi-/3” and “others”). As the topological complexity decreases, 
both Sumax and FI distribute to lower values. The relation found here clearly 
indicates that the flare activity of ARs depends not only on the umbral area 
but also on the magnetic complexity. Moreover, it suggests that the large total 
content of magnetic flux represented by the umbral area is a necessary condition 
and significant magnetic complexity may be another more stringent necessary 
condition for the activation of strong flares. 

3 . 2 . 4 . Twist and Writhe of 6-Spots 

When we classified the 31 ARs by their emergence mode, we found 5 ARs as 
“top-to-top”, 13 as “writhed”, 7 as “quasi-,5” and 6 as “others”. Although the 
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“quasi-/3” and “others” groups hardly ever showed the signs of twist and writhe, 
the “top-to-top” and “writhed” ARs clearly showed the twist and writhe. In 
most cases of the “writhed” type, the signs of twist and writhe agree with each 
other at least in their initial phase (12 cases out of 13). This result strongly 
suggests that the formation of the (5-spot is not due to the gathering of disjointed 
EFRs but due to the emergence of a singly connected structure (left panel in 
Figure EH)) . The “top-to-top” ARs also show high probability of agreement in 
their helicity signs (three cases out of five). We speculate that the “top-to-top” 
type has a more developed writhing structure than the “writhed” type (right 
panel of Figure [20]). However, according to the study of IPark et oi.l (I2012|) . some 
“top-to-top” ARs may be associated with the injection of opposite-sign helicities 
through flux emergence. Therefore the model of “top-to-top” ARs is at present 
controversial and further investigation is necessary to finalize the model. 

The majority of “writhed” ARs were of DK type (12 cases out of 13). The 
UK type was quite rare. The dominance of DK seems natural because a simple 
writhed magnetic tube with upward kink will emerge through the photosphere 
as a normal /3-type magnetic distribution and will rarely be identified as the 
(5-type sunspots. 
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Figure 19. Histograms of umbral area maximum (left) and FI (right) in the 5-state period. 

From the second row downward are shown “top-to-top”, “writhed”, “quasi-/?”, and “others” 
types, respectively, and the top row shows the total. 

3.2.5. 6-Duration and Emergence Type 

As shown in the left panel of Figure [HI more magnetically complex ARs have 
longer i5-duration. We already confirmed that more magnetically complex ARs 
have larger maximum umbral area in Figure (TH] 

In general, larger sunspots have longer lifetime ( jGnevyshevl 1191181 IWaldmeierl 
I1955|l . Recurrent sunspot groups are found to exhibit a slightly larger lifetime 
than that based on the Gnevyshev-Waldmeier relationship by jPetrovay and van Driel-GesztelyH 
(I1997p . It is also well known that larger sunspots tend to be darker and have 
stronger magnetic fields ([Solankil 120021 |Llvingston[l2002p . Furthermore, [Shi and Wan^ 

(I1994P reported that the productivity of (5-spots for X-class flares is closely 
correlated to the lifetime of (5-spots. Our result suggests that the long lifetime of 


SOLA: DELTA_in_Initial_phase_3rd_revised.tex; 24 July 2015; 1:05; p. 25 

















































































































































































Takizawa and Kitai 



Figure 20. Schematic models for ARs. The left panel shows a general model for the 
downward-knotted type such as NOAA 10826. The right panel shows a possible model for the 
“top-to-top” type such as NOAA 10720. 


active 5-spots is closely related to their magnetic complexity (twist and writhe). 
Although magnetic tubes are eventually dispersed due to the action of turbu¬ 
lent convection ( |Petrovay and Moreno-Insertis[ m)97|l and/or fragmentation by 
fluting instabilities, the twisted tubes are more resistant against these dissolving 
actions. It seems difficult that the magnetic lines of force of the tube above the 
solar surface will be torn off by convective motions in lower layers, as the lines 
of force are entangled and hooked to the magnetic tube in the corona. Further, 
the twist will suppress the fluting instability of magnetic tubes (IPriestl r2()14[) . 
Threfore, we argue that the longevity of an active 5-spot depends on its magnetic 
twist. 

3.2.6. Possible Power-law Relations among Parameters 

On every scatter plot (Figures [18] and [21]), we performed the power law fitting 
to the data. Two methods were used to derive the power-law index k. The first 
is the classical least-squares fitting (C-LS) to log-log plots and we obtain the 
power-law index kc. The second is the non-linear least-squares fitting (NL-LS) 
method based on the Gauss-Newton algorithm Icf. IHansen et Qi.1l2()12p . leading 
to /cnl- For the data points in the left panel of Figure [18] (S'umax vs. FI in the 
observed period), we found kc = 1.21 ± 0.24 and A:nl = 1-61 ± 0.33. In the right 
panel of Figure [TH] (^umax vs. FI in the 5-state), we found kc = 1.30 ± 0.31 
and /cnl = 1-87 ± 0.50. FI follows a power-law dependence on 5'umax with an 
index around 2, especially in the 5-state. It is worthwhile to note that [Magara] 
(I2014|) discussed a power-law relationship between the magnetic free energy and 
the total magnetic flux for twisted emerging loops with an index of 2 from his 
numerical simulation. 

We next show the relation between the 5-duration and FI in the 5-state (left 
panel of FigureEI]). Although we could not get the converged solution of NL-LS, 
we obtained kc = 1.54 ±0.34 when we excluded the “others” group as they have 
large dispersion. For the relation between the 5-duration and S'umax in the 5- 
state (right panel of Figure!^, we found kc = 0.76±0.16 and /cnl = 1.07±0.28. 
Combining the last two relations, namely both plots of Figure [H] we can infer 
a power-law relation FI ~ from C-LS which is consistent with the result 

derived in the previous paragraph. 
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Figure 21. Flare Index (left) and maximum umbral area (right) in the (5-state as a function 
of the 5-duration. The symbols are the same as in Figure fTSl Six ARs which showed no flare 
activity in the 5-state were excluded. The error-bars were assigned by considering the effect 
of data gaps (see text). The results of non-linear least-squares fitting (power-law index /cnl) 
and classical least-squares fitting (power-law index kc) are shown by the solid and the dashed 
lines, respectively. 


4. Discussion 

By paying attention to the initial phase of ARs in Cycle 23, we investigated 
what are the important configurations which lead to the flare activity. From our 
investigation, we discovered three representative topological types which have 
strong correlations with the flare activity, namely “quasi-,5”, “writhed”, and “top- 
to-top”. We confirmed that the ARs with higher topological complexity show 
stronger flare activity. 

It is interesting to note that almost all ARs classified as the “writhed” type 
are of DK type. Once the downward knot structure is formed, a dense plasma 
“reservoir” is formed in the knot and the whole structure of the flux tube 
may become stable owing to the effect of the mass as an anchor (left panel of 
FigureEOl). Under a different situation, [Magara and Longcope] (pilOSp proposed a 
similar idea in their numerical simulation of the emergence of a twisted flux tube 
in the convection zone. Subsequently, |Magara| (I2()ll|) has studied the case of the 
emergence of a U-loop formed below the photosphere. |van Driel-Gesztelyi et al] 
(|2000p have reported that the emergence of a U-loop which connected the oppo¬ 
site polarity legs of two U-loops below the photosphere. Recently, IPoisson et al.\ 
(I2013|l studied the possibility of a downward kinked structure in the mid part of 
NOAA10314. Thus the idea of DK structure is not so exceptional and applicable 
to twisted magnetic tubes, with several evidence. The DK structure as an anchor 
may contribute to the stability of the AR. 

Subsequently we checked the case of the “top-to-top” group. They show clear 
twist and writhe structures with common handedness (three cases out of five). 
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Therefore, we infer that the three-dimensional collision is not an accidental one 
but is attributed to the extreme development of a writhing structure. For exam¬ 
ple NOAA10696, classified as the “top-to-top” type, had the second strongest FI 
in our sampled ARs and showed both characters of three-dimensional collision 
and downward knotted structure. The AR is namely a hybrid type of “writhed” 
and “top-to-top”. Such a configuration of this AR reminds us that there may be 
continuity of developments between the “top-to-top” and “writhed” structures. 

We propose a schematic model for NOAA10720 with a highly writhed structure 
as shown in Figure EOl (right panel). We can point out that a similar, highly 
writhed model for NO A A 7926 was proposed bv IPevtsov et aZ.I(ll998p . 

IBrown et al.\ (I2()()d|l studied the rotation of seven sunspots using high resolu¬ 
tion white light images with TRACE. They proposed that a possible mechanism 
for sunspot rotation is the emergence of a pre-twisted flux tube. On the other 
hand. r^n et a7l(|2009p claimed that the rotation of sunspots originates from the 
photospheric flows, based on their study of combined photospheric and coronal 
observations. In our study, the rotation of sunspots themselves and the tilt 
angle rotation of EFR are rarely observed in small and flare-inactive ARs. If 
the rotation is mainly attributed to photospheric flows, rotational motions of 
sunspots will be observed commonly regardless of the sunspot size. However, as 
shown in our statistical study, the ARs of prominent twist and writhe had larger 
area and showed stronger flare activity among the sampled ARs. In addition, as 
described in Section 13.2.51 the pre-twisted flux tube model is in harmony with 
several features of 5-spots. Therefore we conclude that the origin of sunspot 
rotation is the emergence of the pre-twisted flux tube. 

It is known that twist and writhe helicities transform each other in an ideal 
thin flux tube ( |Calugareanu|II95^IMoffat and Riccalll992llRiccallI9951ITorok et al.[ 
I2()l()ll . In this study, we adopted a picture in that the initial twist of ARs will 
be mostly generated by turbulent flows in the convection zone. When the twist 
increases in an AR, the twist transforms into the writhe in the middle part of 
the flux tube by the kink instability (|Tor6k et al.i I2010p . Then it is natural to 
hnd the sign relationship as seen in this study. The previous studies have shown 
that the kink instability facilitates the formation of 5-spots {e.g. ILinton et al.[ 
I1999|l . ITian et ’al] (|2005|1 pointed out that the signs of twist and writhe tend to 
agree in well-developed ARs under kink instability. 

Regarding the classification of 5-spots, we review the relations between the 
types classified by our study and the types by |Zirin and Liggett] (II987|1 as in¬ 
troduced in Section I. We can consider that the hrst type of Zirin and Liggett 
(ZL) is the same as the “top-to-top” type of ours. The third type of ZL may be 
equal to the “writhed” group, although we believe that the configuration is not 
formed by different dipoles but by a connected structure below the photosphere. 
However, the second type of ZL was seldom observed in our study. This type 
was seen in only one or two cases in the “others” group of our study. This type 
may be more likely to occur in recurrent sunspots rather than in sunspots in 
their young stage studied here. The “quasi-/3” type was not considered in the 
ZL classification as they may not produce notable flares. 

We discussed the possible power-law relations among the /375-spot parame¬ 
ters. Our physical interpretation is as follows. FI is a good proxy of the total 
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energy released in flares in an AR and thus is a measure of the total magnetic 
free energy of the region supplied by the gas flow field. The maximum umbral 
area iSumax is a measure of the magnetic flux of the region. The functional 
relation of a power law sometimes is interpreted as the evidence of a scale-free 
process (|Shirnizul 119951 [Shimojo and Shibata 119991 INishizuka et al\ I2009P . In 
our case, the total free energy supply or injection can be considered as due 
to a kind of scale-free processes, probably the helicity injection process by the 
turbulent convective flows or the differential rotation. The process seems to work 
independent of the total amount of magnetic flux. The difference in FI among 
the emergence-type groups may be the difference in topology of the magnetic 
tubes. While the “quasi-/3” type may have a simple form with less twist and 
less free-energy, the “writhed” type may have a singly knotted or kinked tube 
with more free-energy, and the “top-to-top” type may have a doubly or multiply 
knotted tube with much more free-energy. 


5. Conclusions 

We have studied the initial evolution of 31 ARs and derived the following 
conclusions. 

i) From the point of view of topology, emerging (3^5 ARs can be classified broadly 
into three categories: “quasi-/?”, “writhed”, and “top-to-top”. The “top-to-top” 
type has the most complex topology and the “quasi-/?” type is the simplest type. 
The “Writhed” group has medium complexity among the three types. 

ii) ARs of more complex topology tend to exhibit higher flare activity. 

iii) The signs of twist and writhe both tend to be consistent with each other in 
the “writhed” type. 

iv) The downward-knotted structure in the mid portion of the flux tube is the 
essential element of active ARs. 

v) The flare activity of (3^5 ARs are highly correlated not only with the sunspot 
areas but also with the magnetic complexity. 

vi) There is a possible scaling-law between FI and 5'umax- 

Finally, we will give some comments on the future extension of our study. 
Direct measurements of magnetic helicity will be essential to study the free- 
energy accumulation in ARs. Observations of the ARs throughout their lifetime 
along with stereoscopic views from space will be desirable. Further it will be an 
important challenge to diagnose the magnetic topology of ARs in the convection 
zone by improving local helioseismology. 
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